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Abstract 

Incoherent surface scattering yields a statistical description of the surface, due to the ensemble averaging 
over many independently sampled volumes. Depending on the state of the surface and direction of the scat- 
tering vector relative to the surface normal, the height distribution is discrete, continuous, or a combination 
of the two. We present a treatment for the influence of multimodal surface height distributions on Crystal 
Truncation Rod scattering. The effects of a multimodal height distribution are especially evident during in- 
situ monitoring of layer-by-layer thin-film growth via Pulsed Laser Deposition. We model the total height 
distribution as a convolution of discrete and continuous components, resulting in a broadly applicable pa- 
rameterization of surface roughness which can be applied to other scattering probes, such as electrons and 
neutrons. Convolution of such distributions could potentially be applied to interface or chemical scattering. 
Here we find that this analysis describes accurately our experimental studies of (001) SrTi03 annealing and 
homoepitaxial growth. 

PACS numbers: 61.10.-i, 61.10.Kw, 68.37.Ps, 68.47.Gh, 68.55.Ac, 81.15.Fg 

Keywords: Crystal truncation rods, film growth, pulsed deposition, surface roughness, SrTi03, surface scattering, 
X-ray diffraction, X-ray scattering 
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I. INTRODUCTION 



The interaction of a scattering probe with a smooth, crystalline surface gives rise to streaks 
of intensity in reciprocal space known as Crystal Truncation Rods (CTRs). CTR intensity is 
sensitive to atomic-scale surface roughness such as the oscillatory surface roughness that arises 
during the so-called layer-by-layer 1 or polynuclear growth 2-3,4 . Reflection High Energy Elec- 
tron Diffraction (RHEED) is one example of a technique that is often used for in-situ monitoring 
of film growth via CTR intensity . 5 ' 6 i 7 i 8 - 9 Recently, the large monochromatic flux available from 
synchrotron sources (~ 10 13 photons/second) have enabled time-resolved in-situ X-ray scattering 
studies of film growth as well . 10 ' 11 ^ 13 ' 14 ^^ 16 ' 17 ! 18 ^ 9 ' 20 These in-situ scattering techniques enable 
materials engineering at the atomic level via the carefully controlled deposition of sequences of 
single atomic layers . 21122 ' 23124 

Growth of complex oxide thin films with atomic-scale control via scattering-based in-situ mon- 
itoring has recently drawn intense interest. This is due in part to the broad range of materials 
properties that are manifest with changes in stoichiometry, even within a single structure family. 
Pulsed Laser Deposition (PLD) is an attractive growth technique for such epitaxial complex oxide 
films and heterostructures, due to its ability to transfer complex stoichiometrics from the target to 
the film. Many interesting phenomena have recently been reported, including metallic behavior at 
the interface between two band insulators in SrTi03/LaA103 heterostructures 5 , superconductivity 
in BaCu02/SrCu02 heterostructures 25 (neither BaCu02 nor SrCu02 are superconducting in bulk), 
ferroelectricity in SrZr03/SrTi03 heterostructures 26 (neither SrZr03 nor SrTi03 is ferroelectric in 
bulk), and enhanced ferroelectric polarization in CaTi03/SrTi03/BaTi03 heterostructures 27 . 

Despite the increasing popularity of scattering probes as a method for investigating growth, the 
problem of directly comparing models of growth kinetics to scattering data remains a challenge. 
Such a comparison requires descriptions of the surface and the resulting scattered intensity that 
are complex enough to capture the physics of the problem, yet simple enough to lend themselves 
to, for example, least squares fitting. The influence of surface roughness on X-ray intensity can be 
readily calculated for Gaussian roughness 28 or discrete roughness on the order of the out-of-plane 
lattice spacing^ 9 - such as arises during nearly ideal layer-by-layer growth, but we require a more 
generally applicable model to analyze the evolution between these two limits. Previous reports 
have modeled continuous or discrete roughness independently of one another, 30 but real crystal 
surfaces have a miscut, and therefore have concomitant discrete and continuous components of 
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surface roughness. 

Here, we treat the overall surface roughness as a convolution of discrete binomial and contin- 
uous Gaussian distributions, which results in a simple parameterization that is broadly applicable. 
We show that the X-ray intensity along the CTR associated with reciprocal lattice point G can be 
written in a simple, closed-form solution as 



/G ^-^-4Ml-,)shX e±C ' 
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where c' is the discrete step height, n and p are the usual parameters in the binomial distribution, 
and o c is the continuous RMS roughness. We present the application of this model to our studies 
of in-situ X-ray scattering during annealing of (001) SrTi03 and subsequent homoepitaxial growth 
via PLD, and find that our data is accurately described by a simple model for the dependence of 
o c , n and p with film thickness. 

II. THEORY 

In this report, we present a closed form solution of the kinematic scattering theory for surfaces 
with roughness that is continuous, discrete, or a combination of the two. We begin by assuming 
that the first Born approximation is valid, i.e. scattering geometries that do not coincide with 
Bragg conditions. The scattered intensity can be written in the form I(Q) = C(Q) |A(Q) | 2 , where 
C(Q) is a prefactor associated with the scattering geometry, polarization, volume of the unit cell, 
and scattering cross-section . 3 *i 32 i 33 Assuming no surface reconstruction, the scattered intensity for 
a miscut crystal can be calculated by explicitly performing the kinematic sums: 
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where HQ is the momentum transfer, F su b(Q) is the structure factor of the substrate, cq is the out- 
of-plane lattice parameter, S|| is a vector in the plane of the average surface, R(sii) is the position of 
the surface, E{s||} * s a sum over surface positions in each coherently illuminated region, and Y^n ls 
a sum over coherent regions. X-ray scattering measurements can be classified as either coherent, 
partially coherent, or incoherent, depending on the size and coherence volume of the X-ray beam. 
Coherent scattering probes a single realization of the surface (N = 1), while incoherent scattering 
probes the statistical behavior of an ensemble of realizations. This report will focus on incoherent 
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scattering measurements. A simulated reciprocal-space map of the incoherent scattering intensity 
for an (001) -oriented SrTi03 crystal with 0.198° miscut is presented in figure U]-(a). When the 
terrace size is small relative to the transverse coherence length, each reciprocal lattice point G has 
an associated CTR oriented parallel to the surface normal of the average surface. Figure U]-(b) 
compares experimental K-scan data (H = 0, L = 0.5) with the model. The dominant peaks are 
associated with the (000) and (001) CTRs, and satellites can be observed which correspond to the 
(001) and (002) rods. Note also the weak diffuse scattering, which is associated with irregularities 
in the terrace size and thermal diffuse scattering. 

While such brute-force calculations are illuminating, they are quite computationally intensive. 
Alternatively, statistical approaches can be used. Assuming translational invariance, it has previ- 
ously been shown that the scattering intensity of a CTR associated with reciprocal lattice point G 
can in general be written for any surface height distributio n 34 ' 35 

{silHsjj} 

where S|| — Sn is the in-plane vector between scatterers, — = (Q — G) • n, h = r ■ n is the 
surface height as measured along the surface normal (as opposed to the terrace normal), and 

g -»(fii.- G x)[ A (s||)-A(S||)]\ ^ difference function. 



Previous theoretical investigations, primarily for two-level systems such as ideal layer-by-layer 
growth on crystals without miscut, have shown that the total scattering intensity can be split into the 
sum of CTR and diffuse components I tota i = I CTR -\-I di j-j- use 2^^M^l This statistical interpretation 
of scattering has enabled extensive studies of thin-film growth. 30 

Here we focus on the CTR scattering intensity associated with reciprocal lattice point G, which 
may be either specular or off-specular. Given Aft(sn) = fc(sii) — (h), we assume delta-correlated 
height fluctuations (Ah(si\)Ah(s',,)} = o 2 S(sii — Sn) such that 
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-KQ±-Gi_)h{ S]l 



2 



(3) 



For a Gaussian surface height distribution, the scattering intensity can be simplified to yield the 
familiar result 

/G(Q )= /G (Q)e -(ex-G ± )V (4) 

For a flat crystal without miscut, the ideal intensity is 35 



For a crystal with a miscut that yields a surface normal along n, the ideal intensity (which neglects 
the surface roughness associated with the miscut itself), is 



An Atomic Force Microscopy (AFM) image of a (001) SrTi03 surface is shown in figureEKa). 
The existence of a miscut, with irregular terrace sizes, produces a Gaussian height distribution, and 
the resulting X-ray intensity can therefore be described by equation |4j Many surfaces, however, 
have more complicated height distributions. PLD, for example, results in the nucleation, growth, 
and coalescence of many small, two-dimensional islands on the surface. Previous in-situ studies of 
PLD growth have shown that material that deposits on one of these small islands tends to diffuse to 
the step edge and incorporate into the underlying layer , 1 ^ 11 ' 1 ^ enabling extended smooth growth. 
A simulation of the surface shown in figure ERa) after depositing one half monolayer via layer- 
by-layer growth is shown in figure ERb). The system has evolved such that a bimodal distribution 
exists. Therefore, the surface height distribution has both continuous and discrete character. The 
discussion that follows is motivated by the need of a model that can handle the effect of such a 
surface height distribution on the scattered intensity. 

In order to proceed, it is instructive to consider a special case: an ideally flat crystal with zero 
miscut and initial CTR intensity Iq(Q). We then assume random, pulsed deposition where mate- 
rial does not deposit on the same site twice in a single pulse, and no surface diffusion occurs after 
adsorption. The surface height distribution is then given by the well-known binomial distribution: 



k W 

where n is the number of pulses, p is the fraction of the surface covered with new material during 
a single pulse, c' is the step height and k is an index of overlayers. The mean is (h) = npc' and 
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the RMS surface roughness is = c' \Jnp(\ — p). The expectation value in equation |3] can be 
simplified, using the binomial theorem: 




(8) 
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The scattering intensity from a surface with a binomial height distribution is therefore 

/G (Q) l\ wi \ • 2f(Q±-G±W" 
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Alternatively, this result can be derived by performing a sum of the overlayers, as described by 
Robinson, 39 and illustrated in appendix|Aj There it is shown for the more general case where the 
amount of material that adsorbs during each pulse may vary, that the intensity after n pulses is 
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This more flexible description of the surface comes at the expense of additional parameters. In 
figure |3 three CTR intensity profiles for (001) SrTi03 are simulated: an ideally flat substrate; a 
sample with n = 1, p = 0.5, = 1.95A; and a sample with n = 25, p = 0.01 and Cd = 1.94A. 
The CTR intensity is most sensitive to the surface roughness in the so-called anti-Bragg scattering 
geometries, where <2_i_ = v^wi = 1,3,5,... Figure G] illustrates that the RMS surface roughness 
itself is not always sufficient to describe the surface statistically, since the two surfaces with nearly 
identical total roughness have very different CTR profiles. The n,p parameterization provides the 
additional information necessary to model discrete surface roughness on the order of the out-of- 
plane lattice spacing. 

Equation |9] can be rearranged, and by substituting n = o^/c' 2 p(l —p), we arrive at an expres- 
sion that describes how the scattered intensity varies with o^, Q z and p: 
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In the limit of Q± — > G±, equation [TT] can be approximated using first order Taylor expansions, 
sin 2 (jc) ~ x 2 /2, log e (l —x) ~ —x. First expanding sin 2 (x), the expansion of the resulting log 
yields a simplified form of equation fTTl which is independent of p and identical to the Gaussian 
roughness dependence in equation HI 

Alternatively, in the limit of p — ► 0, the log in equation[TT]can be approximated regardless of the 
value of Q z . In this limit, the binomial distribution exhibits Gaussian behavior, and the intensity 
can be written 
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Equation [T2l is applicable to continuous deposition methods, as well as the limit where scattered 
intensity becomes insensitive to the discrete nature of the surface roughness. For experimental 



methods that are insensitive to the crystal miscut, equation [12] therefore also represents the appro- 
priate model for continuous roughness, and is applicable anywhere along the CTR. 

Having derived an expression for the influence of discrete roughness, we return our attention 
to the case where the surface height distribution is a continuous, multimodal function. The con- 
volution of a binomial distribution defined by n = 1 and p = 0.5 with a Gaussian distribution will 
produce the surface height distribution shown in figure El-(b). Using the convolution theorem, we 
now have a solution for the effects of the total surface roughness on the CTR intensity. The rough- 
ness has been parameterized into continuous and discrete contributions, which add in quadrature 
to yield the total surface roughness. For experimental methods that are insensitive to the crystal 
miscut, the total intensity including continuous and discrete contributions of roughness is written: 
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For miscut crystals where the intensity of a single CTR is measured, the total intensity is 
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HI. EXPERIMENTAL RESULTS 



We used surface-sensitive X-ray scattering to study the (001) SrTi03 surface during annealing 
and homoepitaxy at the Cornell High Energy Synchrotron Source. Experiments were performed 
at the G3 hutch on the G-Line beamline. A double crystal synthetic multi-layer monochrometer 
with 1.5% energy bandpass selected 10 KeV X-rays from the 48-pole wiggler spectrum, yielding 
3 x 10 11 photons/s/mm 2 at the sample. The spot size on the sample was set with 2 mm vertical by 
0.5 mm horizontal slit immediately before the sample and the detector resolution was set with a 2 
mm vertical by 0.5 mm horizontal slit at 600 mm from the sample. 



A. Surface Roughening and Annealing 

In order to grow epitaxial films of complex oxides, it is necessary to elevate the substrate tem- 
perature in order to enhance diffusion and reorganization of the adsorbed species. We have found 
that the (001) SrTi03 surface roughens during annealing to 300°C. Figure |4] shows experimental 
specular CTR intensity data from two (001) SrTi03 samples under vacuum, one at room tempera- 
ture (a), the other annealed at 300°C (b). Both samples received a buffered HF etch treatment prior 
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to measurement. The HF etch yields an atomically smooth I1O2 -terminated SrTi03 surface, with 
step edges that are one unit cell in height. During the heating process, we have observed a decline 
of CTR intensity in the (OO5) scattering geometry, which we attribute to an increased step edge 
density that may be related to a loss of oxygen from the SrTi03 surface. 41 A least-squares fit for 
the room temperature data was performed using equation [121 yielding o c = 1.52 A. The CTR 
intensity of the room temperature sample does not appear to be influenced by a discrete surface 
height distribution, as we would expect based on AFM results. The CTR intensity of the 300°C 
data, however, is very strongly influenced by discrete surface roughness; a least-squares fit was 
performed using equation[J21 where n = 1, p = 0.406. The fit of the 300°C data yields a contin- 
uous roughness o c = 0.715 A and a discrete surface roughness of Od = 1.917 A, yielding a total 
surface roughness of a = 2.05 A. 

The temperature-dependent anti-Bragg scattering intensity for (001) SrTi03 samples heated 
in 1 x 10~ 6 Torr O2 through 0.3 Torr O2 is presented in figure As the temperature increased, 
the scattering intensity initially decreased due to increased surface roughness. The temperature 
corresponding to minimum scattering intensity varied with O2 pressure: ~ 250°C for 0.3 Torr, « 
350°C for 1 x 10~ 3 Torr, and « 600°C for 3 x 10~ 6 Torr. As the temperature was increased futher, 
the scattering intensity increased as the surface roughness decreased. Increasing the temperature 
in 0.3 Torr O2 resulted in improved surface roughness at lower temperatures. The baseline error of 
the roughness determination in figure|5]is on the order of 10%, which dominates the determination 
of the errorbars. These results indicate that the surface roughness of (001) SrTi03 can be improved 
by annealing at elevated temperature, and that ramping the substrate temperature in an appropriate 
O2 pressure results in improved surface quality at lower temperatures. 

In figure|6l an (001) SrTi03 substrate was heated to 800°C in 1 x 10~ 6 Torr O2. The scattered 
anti-Bragg intensity was monitored as a function of time, showing modest improvement of the 
surface roughness. After one hour, the O2 pressure was increased to 1 x 10~ 3 Torr, resulting in 
a more rapid improvement of the scattered intensity and surface roughness. The total surface 
roughness in figures |5] and |6] does not approach zero, but a value corresponding to the roughness 
associated with the miscut. This result illustrates that at elevated temperatures, the use of an 
appropriate O2 pressure causes the discrete roughness component to decrease at a faster rate, and 
that both continuous and discrete roughness should be considered during analysis. 
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B. Epitaxial Growth 



Given an acceptable starting surface of (001) SrTi03, we can focus on homoepitaxial growth 
via PLD. A pulsed KrF excimer laser, with 248 nm wavelength and 30 ns pulse duration, was 
focused to provide a 200 MW/cm 2 power density at the surface of a single crystal SrTi03 target. 
The repetition rate was 0.03 Hz, and the growth rate was 0.08 monolayer/pulse. The substrate 
surface temperature was 550°C, and growth proceeded in a 0.02 Torr O2 environment. Time- 
resolved X-ray CTR intensity measurements were made in the anti-Bragg scattering geometry 
in order to maximize sensitivity to surface roughness. PLD growth of (001) SrTiC>3 in these 
conditions proceeds in a layer-by-layer growth mode, where many small two-dimensional islands 
nucleate, grow, and coalesce to form a complete overlayer. This type of growth gives rise to a 
strong discrete component of the surface roughness, which must be accounted for in order to fit 
the resulting oscillatory CTR intensity seen in the experimental data in figure Ufa)- 

The parameterization of roughness contributions in equation [141 enables a direct and simple 
method for analyzing scattering data from in-situ studies of film growth. The binomial distribution 
is used here to describe the state of the surface, not the physics of the adsorption and diffusion 
processes. In place of p, the fractional coverage per pulse, we use to define the fractional 
coverage of the growing layer. For ideal layer-by-layer growth, we set n = <3 2 d /Q{ \ — 0)c /2 = 1 
always, and allow the coverage parameter to vary like a sawtooth waveform, where < < 1. 
Since the scattered intensity depends on the quantity 0(1—0) and not itself, we can equivalently 
define as a triangular waveform 



rico cos(27i/'{fc)— ri) 
L j=0 (2;+l) 2 + 1 
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where < < 0.5 and (k) = (h)/c' is the normalized film thickness. In practice, neither layer 
completion nor half coverage occur at a specific moment across the entire sample, due to small 
position-dependent variations in the deposition rate. The series can be truncated at the first term 
to simulate this effect: = cos(2rc(fc) - %) /4.9352 + 1 /4. 

As the film grows thicker and the continuous roughness increases due to the terraces become 
more and more irregularly shaped, we should expect that the influence of the discrete component 
will diminish. The decaying discrete roughness contribution is modeled using an effective cover- 
age parameter 0' = Y0e~^/ K , where y = 0.8 and K = 18. The presence of a characteristic thickness, 
K, associated with the influence of discrete roughness provides a useful figure of merit for layer-by- 



layer growth. Using scaling arguments described elsewhere, • we have modeled the continuous 
roughness as a power-law o c = a(k)$, with a = 0.617 and |3 = 0.34. These models of the discrete, 
continuous, and total surface roughness as a function of film thickness is presented in figure|7tb). 
By considering both the discrete and continuous components, we obtain good agreement between 
experimentally measured scattered intensity and the model as seen in figure |7t a). 

We have not observed a clear bimodal distribution in AFM data, even for post-deposition mea- 
surements of samples with many small two-dimensional islands or holes on the surface. In figure 
ERb), we show a hypothetical surface that exhibits a bimodal distribution, but we did not account 
for the effect of convoluting the AFM tip with the simulated morphology. Our AFM tips have a 
radius of approximately 10 nm, so it is not surprising that the AFM did not confirm the existence 
of the discrete roughness contribution. In this case, X-rays are more sensitive to such features. 

IV. CONCLUSIONS 

We have used surface-sensitive, incoherent X-ray diffraction for time-resolved studies of (001) 
SrTi03 surface morphology in-situ. Quantitative analysis of specular Crystal Truncation Rod in- 
tensity measurements require a model that is capable of handling the oscillatory surface roughness 
that occurs during layer-by-layer growth, as well as the non-periodic roughness that continues to 
accumulate with increasing film thickness. In general, the RMS surface roughness alone is not 
sufficient to model Crystal Truncation Rod intensity. Crystals with small terrace size relative to 
the coherence length exhibit concomitant discrete and continuous surface roughness, which we 
have modeled as a convolution of continuous and discrete distributions. The total surface rough- 
ness can then be parameterized into continuous and discrete roughness contributions, which add 
in quadrature. We derived a simple, closed-form expression for the scattering intensity which, 
for specific limiting cases, is equivalent to other models existing in the literature. This model is 
capable of handling smooth transitions between roughness that is discrete in nature to roughness 
that is continuous in nature, and therefore relaxes some of the restrictions of existing models of 
reflectivity. 

We have observed good agreement between this model and experimental observations, allow- 
ing quantitative investigations of (001) SrTi03 annealing and thin film growth. We observed a 
surface roughening transition at temperatures that increased with decreasing O2 pressure. The sur- 
face becomes smooth again as the temperature is further increased, and the rate of recovery was 
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improved by increasing the O2 pressure. The maximum surface roughness observed during the 
annealing process was 2.05A, and was primarily discrete in nature. 

The continuous and discrete roughness parameters provide a useful method for analyzing in- 
situ scattering measurements of thin film growth. Growth of oxide thin films via Pulsed Laser 
Deposition tend to proceed with a layer-by-layer mechanism, which gives rise to a strong discrete 
roughness component that is periodic with film thickness. As the film continues to grow, we have 
observed an increase in the continuous roughness as well, which may be due to increasing irregu- 
larity of the terraces or surface relaxation. As the continuous roughness increases, the influence of 
the discrete roughness decreases. The decay rate of the discrete roughness contribution may serve 
as a useful figure of merit for describing layer-by-layer growth. 



APPENDIX A 



We offer an alternate derivation for the effect of a binomial surface height distribution. The 
specularly scattered amplitude from a perfectly flat substrate can be written 

F(Q Z ) 



A ° ~ 1 - e^cn ■ 

The scattering amplitude for an arbitrary surface, just before the n th pulse, is 



A n -i=A +^Q k F(Q z )e- lQ ^ k , 
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where 9^ is the fractional coverage of the W h overlayer. In the case of random deposition, the 
amount of material depositing onto layer ^ — 1 is equal to p(Qk-i — where p is the fraction of 
the total surface covered by a single pulse. The scattering amplitude just after the n th pulse is 



A n = A + £ [p (9n - 0,) + 0,] FiQje-^ 
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= A n —\ -p £ QkF(Q z )e-' &CQk + P t Sk-iF^e-^. 



(A3) 



k=l 



k=l 



The second sum in equation lA3l can be rewritten 
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such that, with some rearrangement, we can solve for /„ in terms of 
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Starting from a perfect surface with no miscut, the intensity after n pulses of random deposition, 
where the fraction p n of the surface covered during each pulse p n may vary, is therefore 
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Figured] (a) Simulated reciprocal space map (H = 0) of intensity for a Ti02 terminated SrTi03 
substrate with 0.198° miscut in the (0K0) direction. The transverse coherence length is assumed 
to be l/jm. (b) Experimental ^-scan data from a (001) SrTi03 substrate, H = 0,L = 0.5. The 
scattering power is normalized to the total incident power. The grey line is simulated with the 
kinematic scattering theory. 

Figure |2j (a) AFM image of a SrTiC>3 substrate. Inset: Histogram of the surface height dis- 
tribution, with a Gaussian fit. The RMS roughness is 1.44 A. (b) Using image (a), simulation of 
random deposition of one half monolayer. 

Figure|3 Calculated intensity along the specular CTR. For n=0, no material has been deposited 
and the substrate is ideally fiat. For one pulse that covers half of the fiat surface with new material, 
complete destructive interference is observed in the anti-Bragg scattering geometries. A nearly 
identical RMS surface roughness is observed for 25 pulses each randomly covering 1% of the 
surface with new material, but the RMS surface roughness alone is not sufficient to model the 
CTR intensity profile. 

Figure |U Reflectivity data for (001) SrTi03 prepared with buffered-HF etch. The orientation 
of the miscut, and the resulting truncation rod splitting, were integrated by the resolution function, 
(a) 25°C, 4 x 10~ 7 Torr. (b) 300°C, 1 x 10~ 5 Torr. 

Figure |3 X-ray scattering studies of SrTi03 surface vs. temperature, (a) (00^) scattering 
intensity vs. temperature, (b) RMS surface roughness vs. temperature, as determined from the 
scattering intensity. 

Figure |6j X-ray scattering measurement of SrTi03 surface vs. time. The sample temperature 
is 800°C, and the O2 pressure is increased from 1 x 10~ 6 Torr to 1 x 10~ 3 Torr after one hour, 
(a) (00^) scattering intensity vs. temperature, (b) RMS surface roughness vs. temperature, as 
determined from the scattering intensity. 

Figure Anti-Bragg X-ray scattering measured in-situ during (001) SrTi03 homoepitaxy via 
PLD. (a) Growth oscillations measured in the (00^) scattering geometry. The model uses the 
roughness depicted in (b), where the total surface roughness includes both continuous and discrete 
contributions, which add in quadrature. 
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